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ON THE CONSTRUCTION OF SUITABLE WEAK 
SOLUTIONS TO THE 3D NAVIER-STOKES EQUATIONS 
IN A BOUNDED DOMAIN BY AN ARTIFICIAL 
COMPRESSIBILITY METHOD 

LUIGI C. BERSELLI AND STEFANO SPIRITO 


Abstract. In this paper we will prove that suitable weak solutions of 
three dimensional Navier-Stokes equations in bounded domain can be 
constructed by a particular type of artificial compressibility approxima¬ 
tion. 


1. Introduction 

Let H C ]R^ be a bounded domain with smooth boundary F = (90 and 
T > 0 be a fixed real number. We consider the three-dimensional Navier- 
Stokes equations with unit viscosity and zero external force: 

dtu — An -|- (n • V) n -b Vp = 0 in x (0, T), 

V-n = 0 innx(0,T), ^ ^ 

and the Navier (slip without friction) boundary conditions for n, namely 

n • n = 0 on r X (0, r), 

(1-2) 

n ■ Dn • r = 0 on F x (0, T), 

where Dn stands for the symmetric part of Vn, n is the unit normal vector 
on F, while r denotes any unit tangential vector on F (Recall also that for 
incompressible fluids An = 2 divDn, to compare with (II.5p in the compress¬ 
ible approximation). 

The system (jl.ll) is coupled with a divergence-free and tangential to the 
boundary initial datum 

n(x,0) = uo{x) on D X {t = 0}. (1-3) 

For the initial value boundary problem dEU-dlSI) it is well-known that, for 
any tangential and divergence-free vector field no S L^(D), there exists at 
least a global weak solution in the sense of Leray-Hopf, see [DElElIiniEniEs] 
for the generalization of the classical results of Hopf [23], which was obtained 
for the Dirichlet case. 
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Remark 1.1. In the sequel it will be enough to assume that F G and, 
to avoid technical complications, we also assume that the domain is simply 
connected and that it cannot be generated by revolution around a given axis 
(this could be relevant in the steady case, when dealing with the symmetric 
deformation tensor and Korn type inequalities), see [5]. 


On the other hand, the uniqueness and the regularity of the weak solu¬ 
tions represents a problem still open and very far to be understood. The 
best regularity result which is available for weak solutions of the three- 
dimensional Navier-Stokes equations is the Caffarelli-Kohn-Nirenberg theo¬ 
rem which asserts that the velocity is smooth out of a set of parabolic 
Hausdorff dimension zero. For example it implies that there are not space- 
time curves of singularity in the velocity. However, the Caffarelli-Kohn- 
Nirenberg theorem holds only for a particular subclass of weak solutions, 
called in literature suitable weak solutions,'' see Section [2] for the precise 
definitions. Roughly speaking a suitable weak solution is a particular Leray- 
Hopf weak solution which, in addition to the global energy inequality, satis- 
hes in the sense of distributions also the following entropy-type inequality: 


dt 


u 


+ V- 


1 , 


u\^ +p] u] — A 


u 


+ |Vup<0. (1.4) 


Inequality ()1.4p is often called in literature generalized energy inequality and 
is the main tool to prove the partial regularity theorems HU [25]. Since 
uniqueness of Leray-Hopf weak solutions (within the same class of solu¬ 
tions) is not known, each method used to prove existence of weak solutions 
may possibly lead to a different weak solution. This in turn implies that it 
is a very interesting problem to understand which ones of the different ap¬ 
proximation methods (in particular those important in applications and in 
the construction of numerical solutions), provide existence of suitable weak 
solutions. This question has been considered for many approximation meth¬ 
ods laEnaiei, and it is worth to point out that the solutions constructed 
by the Leray method |26] and the Leray-a variant turn out to be suitable. 
Particularly important are the results obtained by Guermond [201 |2T] where 
it was proved that some special Galerkin methods lead to suitable weak 
solutions. However, understanding whether solutions obtained by general 
Galerkin methods (especially those ones based on Fourier series expansion) 
are suitable has been posed in [3] and it is still not completely solved. 

In the spirit of understanding whether certain methods used in the nu¬ 
merical approximation produce suitable weak solutions (problem empha¬ 
sized in [22|) in this note we prove that weak solutions obtained as a limit 
of a particular artificial compressibility method are so. More precisely, the 
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approximation system we consider is the following: 

dtu^ + • V) div 'u'" - 2 div Dtt^ + Vp*" = 0 

—e + div = 0 

is coupled with Navier boundary conditions 

■ n = 0 on r X (0,T), 

n • • T = 0 on r X (0, T), 

and a tangential to the boundary initial datum 

M^(x,0) = Uq{x) on O X {t = 0}. 

Concerning the pressure, it is natural to impose on it Neumann boundary 
conditions, 

Qrf 

^ = 0 onrx(0,r), (1.8) 

and normalization by zero average over the domain: 

/ p^{x,t) dx = 0 a.e. t S (0, T). (1-9) 

Jn 

Remark 1.2. The nonlinear terms in (HSI) can be also be written in the 
following equivalent divergence-type form 

nl{u^,u^) := {u^ ■V)u^ + divu^ = div(M^ (g) u^) “ ^ divu^, (1.10) 

where as usual (a® b)ij := aibj, for i,j = 1,2,3. The factor ^u^divu^, 
which vanishes in the incompressible limit, is needed to keep the usual energy 
estimates. 

We recall that artificial compressibility methods were first studied by 
Chorin [nmi] and Temam [31j and are relevant in numerical analysis since 
they relax the divergence-free constraint, which has a high computational 
cost. See also the recent results in [28] for the Cauchy problem with —e Ap^ 
replaced by ep^ (as in [5]) in equations (II.5p . In addition, the case where 
the term —e Ap^ in equations (II. is replaced by e dtp‘" has been considered 
by different authors, see |32l[l5|. Moreover, the convergence when e —)• 0 to 
a suitable weak solution in the whole space was considered in m for the 
scheme with the time derivative of the pressure. 

Here, we focus especially on the fact that we have a domain with solid 
boundaries and the coupling with the Navier boundary conditions makes 
possible to obtain appropriate estimates on the pressure. The main theorem 
of this note is the following. See Section [2] for the notations and the main 
definitions. 


in X (0,T), 

in X (0, T), 

(1.5) 

( 1 . 6 ) 

(1.7) 
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Theorem 1.3. Let Ll be bounded and of class C^’^. Let {{u^be a 
sequence of weak solutions of the initial value boundary problem (iLSD-dLH]) 
with 

div Uq = 0 in Ll, 

and {«o}e bounded uniformly in H^{Q) such that 

Uq —>• uo strongly in L‘^{LI). 

Then, the following e-independent estimate on the pressure holds true 

3C>0: Wp^W , <C, Ve > 0. (1.11) 

"L3((0,T)xO) — ’ ^ ^ 

Moreover, there exists u, a Leray-Hopf weak solution of (fTTTD-ra. and an 
associated pressure p G L3((0,T) x Ll) such that 

weakly* in L°°{0,T;L^{n)), (1.12) 

weakly in L°°{0,T-, L‘^{n)), (1.13) 

weakly in L3((0,T) x 11). (1-14) 

Finally, {u,p) is a suitable weak solution to the Navier-Stokes equations (jl.ip - 

([Llll-dliD- 

The existence of weak solutions of (ll.5l) - (jl.8p is quite standard and can be 
easily proved by smoothing suitably the non linear terms. The main obsta¬ 
cle to prove Theorem 11.31 is to get e-independent estimates on the pressure. 
Since we do not have the divergence-free constraint we cannot deduce inde¬ 
pendent estimates on the pressure by using the classical elliptic equations 
associated to the pressure. Moreover, also methods based on the semigroup 
theory as in |29j seem not directly working here, since the approximation 
system does not immediately ht in that abstract framework. We will get the 
necessary a priori estimates on the pressure from the momentum equations 
as in the case of compressible Navier-Stokes equations [HI [27]. This ap¬ 
proach has been introduced in uni to study a class of non-Newtonian fluids 
and used also in [12] to address the analysis of models for the evolution of 
the turbulent kinetic energy, but without considering the question of the 
local energy inequality. In particular, we point out that it is in the pressure 
estimate that we need to employ the Navier boundary conditions, since they 
allow us to control the various term arising in the integration by parts. 


£ * 

u ^ u 
Vtt^ ^ Vu 

P^ ^P 


Plan of the paper. The plan of the paper is the following. In Section [2] 
we recall the main definitions and the main tools we will use in the proof of 
Theorem [L3l In Section [3] we prove the main a priori estimates needed in 
the proof of Theorem 11.31 and finally in Section |4] we prove Theorem 11.31 
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2. Preliminaries 


In this section we fix the notations and we recall the main definitions and 
tools we will use to prove Theorem 11.31 Given C the space of C°° 
functions or vector fields on If tangential to the boundary will be denoted 
by We will denote with the standard Lebesgue spaces and 

with II • lip their norm. The classical Sobolev space is denoted by 
and its norm by || • ||fc^p and when k = 1 and p = 2 we denote with 

As usual we denote the L^-scalar product by (., .). Finally, we 
denote by and the space of the tangential to the boundary 

and divergence-free vector fields respectively in L^(n) and By 

we denote the scalar product between vectors, while by “: ” we denote the 
complete contraction of second order tensors. Finally, given a Banach space 
we denote by C^dO, T]]X) the space of continuous function from the interval 
[0, T] to the space X endowed with the weak topology. 

Let us start by giving the precise definition of Leray-Hopf weak solution 
for the initial value boundary problem (jl.ip - (ll.2p . 


Definition 2.1. A vector field u is a Leray-Hopf weak solution to the Navier- 
Stokes equations (irili-D-dol) if the following proprieties hold true: 

1) The velocity u satisfies u G 6*^(0, T ; n T^(0, T ; Ff^,.(n)); 

2) The velocity u satisfies the following integral identity 

j J (^u ■ dtfi — Vu :V(j) — {u ■V)u ■ dxds— Ur u ■ Vn • (pdSds 

= — uq ■ 4>{0) dx, 

Jn 

for all smooth vector fields (f>, divergence-free, tangential to the boundary, 
and such that 4>{T) = 0; 

3) The velocity u satisfies the global energy inequality 


\\\u{t)\\l + \\Vu{s)\\lds- 


u-Vn ■ udSds <-\\uo\\ 2 - (2T) 


f u-'Vn-udS:= f UiUj---^dS. 

Jr Jr (Jxj 


Remark 2.2. In the above definition and in the sequel, we recall that 

3 
u 

T JT ^ '-'•^3 

and we also recall that the boundary condition ()1.2[) are encoded in the weak 
formulation, see also |24j . 

The next definition we recall is that of suitable weak solutions. 


Definition 2.3. A pair {u,p) is a suitable weak solution to the Navier- 
Stokes equations (irili-D-dol) if the following properties hold true: 

1) u is a Leray-Hopf weak solution with associated a pressure p such that 

pe lI{{ 0,T) xQ); (2.2) 
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2} The generalized energy inequality holds true 

/ / \Vu\‘^ (j) dxds < / / —— [(pt + A(p^ dxds 

Jo Jq Jo Jq ^ ^2 3 ^ 

( I I ^ \ 

^- -Vcpdxds, 

for all non-negative (p € C'“((0,T) x it). 

We want to point out that the real difference between the Leray-Hopf 
weak solutions and the suitable weak solutions relies in the local energy 
inequality. Indeed, it is always possible to associate to a Leray-Hopf weak 
solutions u a pressure p (modulo arbitrary functions of time) which belongs 
to the space L3((0,T) x H). This is standard in a setting without physical 
boundaries by using the Stokes system. In the case of Dirichlet boundary 
conditions this was first proved by Sohr and von Wahl in |29] . Finally, in 
the case of Navier boundary conditions this can be deduced again by the 
elliptic equations associated to the pressure p, see [8]. On the other hand, 
it is not known how the deduce the generalized energy inequality (II.4h since 
the regularity of a Leray-Hopf weak solution is not enough to justify the 
chain rule in the time derivative term and in the non linear term. 

We pass now to the precise definitions concerning the compressible ap¬ 
proximation we will consider. 


Definition 2.4. The couple {u^,p^) is a weak solution to the compressible 
approximation (II3])-(IIll)-(in|-(IISl)-(II31) if: 

1) The velocity and pressure {u^,p^) satisfy 

G L°°(0,r;L2(Ll))nL2(0,r;//i(H)), 

G L2((0,T) X H); 


2) The velocity and pressure {u^,p^) satisfy the following integral identity 


fL 


— ■ dt<p + nl{u^, u^) ■ (p + 2Du^ : Dtp -|- Vp^ ■ (p 


-\-£ Vp^ • VV’ + div Ip 


I 

Jq 


dxds = / ^0 ■ dx, 


for all smooth vector fields cp tangential to the boundary, such that <p{T) = 0 
and for all smooth scalar fields ip with zero mean value; 

3) The velocity and pressure {u^,p^) satisfy the global energy inequality 


-\\u^ 

2 " 


-k 


Jo 


-k e ||Vt 


, e ||2 


ds < -lln! 


s\\2 


Vf>0; (2.4) 
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4) The velocity and pressure satisfy the local energy inequality 

pT p r I ^1^ 

/ / (2\Du^\^ + e\Vp^\^\(f)dxds < / / —-— ((/>f + Ac/)) dxds 

Jo Jn^ ' Jo Jn 2 


+ 


/ / p^u^ — ep’^Vp’^ + divu^^ • V(f)dxds (2.5) 

0 t/ 


+ 


/ j u®u : VJidxds, 

Jo Jn 


for all non-negative (j) G C^{{0,T) x 12). 


Remark 2.5. The global and local energy estimates can be rewritten in a 
more useful form by performing some integration by parts. Observe that for 
any v G it holds (cf. [5l Sec. 2]) that 

2 / Du : Dudx = ||Vu|p + II divu|p — / v-Vn-vdS. 

Jn Jv 

Moreover, by a standard compactness argument (see [5l Lem. 2.3]), it follows 
that there exists c = c(D) such that 

2 [ Dv:Bvdx> c||Vuf Vu G 2/^(D), 

Jn 

and the energy inequality (|2.4I) can be rewritten either as follows 

1 ft ^ 

/ llVu^f+ ||divu^f+e||V/"2 

20 

— u'^ ■ Vn • dS 




or even as 


L 

^\\u’{t)f + l‘[4Vu‘f + 4Vp‘ 


( 2 . 6 ) 


ds < -||unlP) 
- 2" oil ’ 


ds < 


(2.7) 


Concerning the local energy inequality (12.51) . with several integration by 
parts we get the following equivalent formulation (which is more similar to 
the standard one for the incompressible case), 

pT p r'^ r I ^1^ 

/ / f|Vu^|^ + I divu'^p + e|Vp^p'))/)(ix(2s < / / ^ {(pt + Ac/)) dxds 

Jo Jn^ 2 Jo Jn 2 


+ 


/ / — ep^Vp^ — u^ divu^^ • Vfdxds, 

t/ 0 t/ 


for all non-negative (j) £ ^“((0,7") x 12). 


( 2 . 8 ) 


The results in m, when specialized to the case r = 2 imply the following 
theorem and the local energy inequality follows from the improved regularity 
of the pressure G L^((0, T) x 12), which is valid (but not uniform) for 
all e > 0. 










BERSELLI AND SPIRITO 


Theorem 2.6. Let be given Uq G L‘^{LI), then there exists a weak solution 
to the compressible approximation CZ])-(lOI)-(1121)-(IIJ)-(lL9l) . 

Finally, we recall the classical Aubin-Lions Lemma which will be useful 
to obtain compactness in time in the passage to the limit from (|1.5I) to (11.11) . 

Lemma 2.7. Let X, B, and Y he reflexive Banach spaces. For e > 0 let 
be a family of functions uniformly hounded in LJ’{0, T ; X) with p > 1 
and let {dtu^}e be uniformly bounded in L^{0,T;Y). with r > 1. If X is 
compactly embedded in B and B is continuously embedded in Y, then {u^}e 
is relatively compact in L^(0,r;i?). 


3. Some estimates on the pressure independent of e 

In this section we prove the e-independent estimate (II.lip for the pressure, 
which represents the most relevant technical part to show the convergence 
towards a suitable weak solution. 


Lemma 3.1. Let {u‘^,p‘^) be a weak solution of system (ll.5l) - (ll.7p - (|1.8p - 
dini). Then, there exists a constant C > 0, independent of s, such that 

sup e||/(t)|||+ [ ||/(s)||| ds < C. 

0<t<T 3 Jo 3 

Proof. Let ce = | and let g'^ be the unique solution (normalized by a vanish¬ 
ing mean value) of the following Poisson problem with Neumann boundary 
conditions 

Ag^ = \p^\°‘~‘^p^ — f \p^\°'~'^p^ dx in n X {t}, 

< (3.1) 

^ = 0 on r X {t}. 

K (jTl 

Remark 3.2. The number a — 2 is negative but the expression we write is 
legitimate since we are not really dividing by zero in sets of positive measure. 
In fact, by using the regularity of p^ from Dehnition 12.41 it turns out that 
the function \p^\'^~‘^p^ is well-defined and belongs at least to L^((0,T) x 12), 
even if not uniformly in e > 0. 


We use now the vector field as test function in the hrst equation of 
the system (11.511 . Note that, Vg'^ is tangential to the boundary and for each 
hxed e > 0 and we have that Xg^ is smooth enough to make the integral 
below well-defined. Let t G (0, T), then we get (recall (ll.lOp i 



-|- nl{u^, u^) — 2 divDu^ -|- Vp^) • Vg^ dxds = 0. 


Integrating by parts the term involving the deformation gradient by using 
the Navier conditions ()1.7|) and the one involving the pressure and by using 








SUITABLE SOLUTIONS BY ARTIFICIAL COMPRESSIBILITY 


9 


the definition of , the Neumann boundary conditions, and the fact that 
has zero average we get 


[ [ \p^\^ dxds = f [ [dtu^ + nl{u^,u^)) 
Jo Jn Jo Jn L 


2„e 


Vg + 2 Dm : V g 


dxds. 

(3.2) 


Let us consider the term with the time derivative. We can use the Helmholtz 
decomposition to write 


= Pu^ + V(?. 


where P is the Leray projector. Note that, since • n = 0 on F x (0, T) it 
holds that n • Vg = 0 on F x (0,T) as well. By using the second equation 
of (jl.5[) we have that = jdiv{Pu^ + Vg) = ^divVg = ^Ag. This 
means that for a.e. t G (0,T) we have in a weak sense 


A = 0 in H X {t}, 

A(p'-|)=o in rx{t}. 

Then, by the uniqueness of the Neumann problem under the usual normal¬ 
ization of zero mean value, we get that 


/ = -• 
e 

Then, we observe that since • n = 0 on F and div Pu^ = 0 in D, then 

n dtu^ ■ Vg^ dxds = — f f {div dfU^) g^ dxds = — f f Adtqg^dxds 
! Jo Jn Jo Jn 



= —e / / Adtp^ g^ dxds. 

lo Jn 


We then integrate by parts twice and by observing that n ■ Vp^ = 0, (hence 
a fortiori also n • Vdtp^ = 0) and n ■'S/g^ = 0 we obtain 


fl 

Jo Jn 


dtu^ ■ Vg^ dxds = —s f f dtp^ Ag^ dxds. 

Jo Jn 


'0 Jn 

Hence, by recalling the definition of g^ via the boundary value problem (13.ip 
we get 

n dtu^ ■ Vg^ dxds =—e / / dtp^{x,s) p^{x,s)\p^{x,s)\°‘~^ dxds 

! Jo Jn 


where we used again the fact that p'^ has zero average for any t G (0, T). 

Then, by recalling that a = |, by integrating in time the last term, and 
taking into account the fact that all the initial data {mq}£ are divergence-free 
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(hence p^(0) = 0) we get that the equation (|3.2p reads now 

^ll/(^)III+ / ll/(s)III = [ [ nl{u^,u^) -Vg^ dxds 
3 3 Jo 3 Jo Jn 

+ 2 / / Du^ : dxds 

Jo Jn 

=: h + h- 

We estimate the first integral from the right-hand side as follows: 


h < 


[ ||re^(^l^u^)||l5 ||V/||i5ds 

Jo 


< [ \\nl{u^,u‘^)\\w\\D‘^g‘^\\5 ds 

Jo 14 2 

< [ ||n/(«^,tt^)|| 15 IIp^III ds 

Jo 14 3 

rt 5 1 /■* 5 

/ 11/Hick'S, 

/o 3 


< C 


f \\nl{u^,u^)\\l5ds + ^ 
Jo 14 ° 


where we have used Holder, Sobolev, and Young inequalities. Finally, by 
convex interpolation we have that if G L°°(0, T; L^(f2))nL^(0, T; 
then the nonlinear term nl{u^,u^) belongs to Ls (0, T; Lw (n)). 

Hence, we have proved that there exists a constant Ci, depending only on 
the L^-norm of the initial data (through the energy inequality (j2.7p ) such 
that 


Ii<C + 




3 


Concerning the term I 2 we have, by using Holder inequality and again the 
elliptic estimates 

rt 


[ [ |Vn^ 
Jo Jn 


72„e| 



5 5 ds 

Jo 

3 2 

/■* 

2 

<c{n) / 1 

Vrt^ 2 / 5 ds 

Jo 

( 

r / 

<c(o,r)^ 


1 f 

t 5 


Wp^Wlds, 

3 


ds + 


f \\P‘ 

Jo 


II ids 


for a constant C 2 , depending only on fl, T, and the L^-norm of the initial 
data (again through the energy inequality (I2.7P '). Collecting all estimates, 
we finally get 

^ii/wiii+^ ri/(s)iitds<ci+c2, 

o 3 ^ Jo 3 
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ending the proof. By the way we also proved that 

3 5 

is bounded uniformly in L°°(0, T; Ls (0)), 
even if this information will be not used in the sequel. □ 


4. Proof of the Main Theorem 


By using the estimates from the previous section we can now prove Theo¬ 
rem oi in an elementary way by using standard weak compactness method. 
By recalling (j2.6p - (ll.llh we have that, up to sub-sequences, there exist 
u S T°°(0,r;L2(j^))nL2(o,r;//i(f2)) andpG Lf((0,r) x Q) such that 

^ Vu weakly in L^{0,T x Q), 

^ u weakly* in L°°(0, T; L^po)), (4.1) 

^ p weakly in La ((0, T) x Q). 

Moreover, by (12.61) if follows that y/s\7p^ is uniformly bounded in T^((0, T) x 
fl), and then we have that £ 2 +^ Vp'" converges strongly to zero for all positive 
6. Hence, in particular, we have 

^ 0 strongly in L'^{{0,T) x fl). (4.2) 


It also follows that u is divergence-free. Indeed, by the second equation 
of (jl.5p we get 

/ / di\ il;dxds = ^/e / < >H-i ds 

Jo Jn Jo ’ ° 


=-v/e f [ y/eVp^ ■Vil}dxds 
Jo Jn 

<^/e[ ||v^Vp^||2||VV'||2 ds 
Jo 

< ■\/c||-\/cVp^||i2(o,T;L2(H))l|VV^IlL2(0,r;L2(f7)). 

By taking the supremum over the functions i/; G L‘^{0,T] Hq{Q,)) and by 
using the estimate on p^ from (j2.6|) we get that 

divM^^O weakly in L2(0,T; ff~^(n)). 


Then, by the uniqueness of weak limits we get that u is divergence-free. 

The next step is to prove the strong convergence of in L‘^{(0,T) x Cl). 
We can now use Aubin-Lions, Lemma 12.71 provided that we can show some 
estimates on the time derivative of the velocity, and these are usual obtained 
by comparison. We observe that p^ is uniformly bounded in Ls ((0, T) x fl)), 

4 

the nonlinear term nl{u^, u^) is uniformly bounded in La (0, T ; H~'^{Ct)), and 
divDu^ is uniformly bounded in L?‘{d,T]H~^{Cl))., hence we get that 

dtvJ G L3(0, T; IT“^’t(n)), uniformly with respect to e. 
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Then, by applying Lemma [2.71 with X = B = L^(0), and Y = 

we get 

^ u strongly in L‘^{{0,T) x 0). (4.3) 

Then, we investigate the convergence of Pu^ and Qu^, where P is the Leray 
projector and Q := I — P. By applying P to the first equation of (ll.Sp we 
get that 

dtPu'^ - P{Au^) + P{{u^ ■ V) u^) - ^ Q divu^)^ = 0. 

From this equation we show, again by comparison, that dtPu^ is bounded in 

4 

L3(0,T; uniformly with respect to e, and since Pu^ is bounded 

in L‘^{0,T; P[l ^{^)) we can apply again Lemma 12.71 to obtain in a very 
standard way that 

Pu^ —)• Pu = u strongly in L^((0, T) x O)). 

Then, we have that 

= 11?^^ — u + u — T’u^||2 

< \\u^ — n ||2 + \\Pu — T’u ^||2 —>•0, as e ^ 0. 

This in turn implies, 

Qu^ —^ 0 strongly in L^((0,T) x 14)). 

The next step is to prove that Qu^ converge to 0 strongly in La ((0, T) x n)). 
By using an interpolation inequality and a Sobolev embedding theorem 
together with the Poincare inequality (valid for tangential vector fields, 
see [18]) we get 

WQu^Wl < \\Quy\l\\Quy\l <C\\Qu^\\l\\XQuy\l 
2 

and due to the definition of — Pu^ = Xq. By integrating in time 

we get 

/ WQu’' is)\\l ds < c [ ||Qu=(s)|||||Vu=(s)||| ds 
Jo 2 Jq 

rp ^ rp 3 

< sup \mml ([ iiQu^(s)iiid.y f [ iivu^(.)iiidsy. 

0<t<T \Jo J \h J 

Then, by using the regularity of {u^,p^) and construction of the Helmholtz 
decomposition, it holds that G L°°(0, T; L^(H)) and we get 

Qu^ 0 strongly in L2((0, T) x 14)). (4.4) 

The last convergence we need to prove concerns the term e Xp^. By consid¬ 
ering the second equation of (|1.5p we get that 

e Ap^ = div = div Qu^. 
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By standard using estimates on the elliptic equations with Neumann bound¬ 
ary conditions, see m, we get the following estimate: 

e llVp'^ll 5 < 

By integrating in time and using (14.4p we get then 

o 


[ llv/III 

Jo 2 


0 . 


By using (I4.ip . (I4.2p . (I4.3p . and (I2.6p it is straightforward to prove that u is 
a Leray-Hopf weak solution. In order to prove that {u,p) is a suitable weak 
solution of the Navier-Stokes equations dEUl-dLl we only have to show is 
that {u,p) satisfies the local energy inequality (12.3p . Since (p > 0, from (j2.8p 
we have that {u^,p^) satisfies 

e|2 


f f \Vu^f(j)dxds < f [ + 

Jo Jn Jo Jn ^ J 

£|2 


+ ( 


u 


u 


+ p^u^ — p^Vp^ — div • Vp 


By weak lower semicontinuity of the L^-norm, the fact that Vn' 
weakly in L^((0,T) x Q), and since p > 0 we have that 

f [ \Vpdxds < ImiiTd f f pdxds. 

Jo Jn ^^0 Jq Jq 

Since ^ u strongly in L^((0,T) x fl), we get 
fT f U£|2 

/ / —{Pt + A0) dxds 

Jo Jn 4 


dxds. 

(4.5) 
^ Vn 


10 Jn 


■ (^pt -|- Ap) dxds, as e ^ 0. 


Next, by interpolation we also have that ^ u strongly in L^(0, T; L^(n)) 
and that is bounded in L‘^{0,T; L^{Q)). Consequently, it also follows that 


fl" 


u 


e\2 


f IttP 

Vpdxds^ / / u——-pdxds, 

Jo Jn 2 


as e 


0 . 


Now, we estimate the last two terms in (14.5p . We start by estimating the 
term involving the pressure. We have that 

p^ ^ p weakly in La (fl x (0, T)), 

while by standard interpolation argument 

^ u strongly in L^{Q. x (0,T)). 

These in turn imply that 


fl 


p^u^ ■ Vpdxds 


/ / pu-Vpdxds, 

Jo Jn 


as e 


0 . 
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Concerning the integral 


A := e 




■ V4>dxds, 


(quadratic in the pressure) we argue as follows: By using Holder inequality 
we get that 

, fT ^\l / rT 

1^1 < Ce 


< Ce 


[ l|Vp' 

Jo 

f l|Vp‘ 

Jo 


,e|| 2 

2 


/ ll/llf 

'0 3 , 


2 

5\ 5 


,e|| 2 

2 


and also A goes to 0 as e —)• 0. 

Finally, the we consider the last term, namely 

B :=— / div ■ ^(pdxds, 

Jo Jn 


and since divn'" converge weakly to 0 in L^(0, T; and is uni¬ 
formly bounded in L^(0, T ; (fl)) we get, by uniqueness of the weak limits, 

that 

divu^ ^ 0 weakly in L^((0, T) x bl). (4.6) 

By using the fact the —)• tt strongly in F^((0, T) x H) we have that B goes 

to 0 when e vanishes, hence that (12.3p is satisfied. 
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